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We discuss designer Hamiltonians — lattice models tailored to be free from sign problems ("de- 
signed") when simulated with quantum Monte Carlo methods but which still host complex many- 
body states and quantum phase transitions of interest in condensed matter physics. We focus on 
quantum spin systems in which competing interactions lead to non-magnetic ground states. These 
states and the associated quantum phase transitions can be studied in great detail, enabling direct 
access to universal properties and connections with low-energy effective quantum field theories. As 
specific examples, we discuss the transition from a Neel antiferromagnet to either a uniform quantum 
paramagnet or a spontaneously symmetry-broken valence-bond solid in SU(2) and SU(iV) invariant 
spin models. We also discuss anisotropic (XXZ) systems harboring topological Z2 spin liquids and 
the XY* transition. We briefly review recent progress on quantum Monte Carlo algorithms, including 
ground state projection in the valence-bond basis and direct computation of the Renyi variants of 
the entanglement entropy. 
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I. INTRODUCTION 

Understanding ground states of interacting quantum 
systems is one of the central defining challenges of con- 
densed matter physics. On the one hand, quantum field 
theory is a powerful and general framework for studies of 
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low-energy properties of complex, strongly correlated and 
entangled quantum matter On the other hand, this 
continuum approach has potential pitfalls due to its gen- 
erality, and limitations due to its insensitivity to lattice 
scale physics, making it important to also carry out direct 
studies of microscopic lattice Hamiltonians. A fruitful di- 
rection is to combine the advantages of both methods by 
analyzing data from unbiased numerical investigations of 
lattice Hamiltonians with quantum held theoretic pre- 
dictions of emergent low energy properties. Such con- 
nections between the two approaches can have signifi- 
cant synergistic effects in advancing our understanding 
of quantum many-body phenomena — as they have had 
and continue to have in classical statistical mechanics, 
especially in studies of thermal phase transitions [H, 0] ■ 

Computational many-body research faces important 
challenges since the quantum models of interest are of- 
ten beyond the reach of existing algorithms. Two- 
dimensional (2D) and 3D frustrated quantum spin sys- 
tems and fermion systems are the main groups of diffi- 
cult systems, since they suffer from the infamous "sign 
problem" in quantum Monte Carlo (QMC) simulations 
|4]-|6(. The sign problem arises when the weight func- 
tion in the configuration space constructed using the 
Euclidean path integral or other mappings to an effec- 
tive statistical-mechanics problem is not positive definite, 
hence invalidating the interpretation of the weights as a 
probability distribution for importance sampling. The 
approach we advocate and review here is to construct 
particular "designer Hamiltonians" , which host interest- 
ing ground states and quantum phase transitions but still 
are amenable to large-scale QMC studies without sign 
problems (i.e., they are designed and "de-signed" to be 
practically useful). 

A common misconception is that interesting designer 
Hamiltonians should not exist — if a model is sign- 
problem free it must be trivial or uninteresting. In this 
article, we will review some recent counter-examples to 
this pessimistic view, from the area of quantum mag- 
netism. Moreover, advances in QMC methods |7H12j| 
have made it possible to truly approach the low-energy 
limit of many important models without any approxi- 
mations (beyond controllable statistical errors), and to 
compute quantities of current interest in condensed mat- 
ter physics and quantum information theory, e.g., the 
fidelity 13 1, the geometric tensor [l4|, and the entan- 
glement entropy [ljl [l6j]. We will briefly address some 
of these technical developments as well. We first make 
some further remarks on field theory and the philosophy 
of the synergistic use of designer Hamiltonians, and then 
outline the topics covered in this Review. 



Field theory and numerical studies of 
Hamiltonians 



The natural starting point for theoretical models of 
condensed matter systems arc microscopic Hamiltonians. 



An effective field theoretic description of a given Hamil- 
tonian is constructed based on a combination of insights 
and assumptions on the properties of the Hamiltonian 
that control the long distance physics according to the 
renormalization group (RG) approach. These properties 
most famously include symmetry and dimensionality, but 
more subtle effects like the nature of the topological de- 
fects also play an important role. The fixed points of the 
field theory control the long distance physics. The RG 
approach may result in flows to a previously known fixed 
point, or to some fixed point whose properties are not 
known. The RG itself often relies on approximations that 
are not easy to justify, and, if the flow is to some unknown 
state, it is challenging to characterize it. Often to control 
the RG flows, analytic expansions around some special 
simplifying limits are considered. Well known examples 
are the e expansion around the upper critical dimension 
and the 1 /N expansion in the number of components iV 
(spin components, flavors of fermions, etc.) [13, EH- The 
extrapolation to the e or N of interest is often problem- 
atic, however, because only low-order expansions are fea- 
sible in practice — or worse still, the expansion may not 
converge at the value of e or N of primary interest. 

One example of the above approach of deriving effec- 
tive field theories from microscopic Hamiltonians by a 
limiting expansion is the (d + l)-dimensional nonlinear 
(T-model description of interacting quantum spins in d 
dimensions, which was derived by Haldane in the semi- 
classical limit of large spin S [l9l — [211 ] . Another example — 
central to this review article — is the non-compact CP^ -1 
description of the intriguing Neel to valence-bond solid 
(VBS) "deconfined" quantum-critical point in 2D S\J(N) 
quantum magnets, which can be justified in the large- N 
limit [H, [23| . Whether or not the theory remains cor- 
rect down to small N, in particular N — 2, remains a 
challenging question. 

In view of the uncertainties with analytic methods 
highlighted above, it is important to test the predictions 
of field theories in some unbiased way, starting from a 
microscopic description of the system or phenomenon of 
interest. In some cases there are exact solutions, e.g., the 
Bethe Ansatz solution of the S = 1/2 Heisenberg chain 
[24j has been crucial for testing the non-linear tr-model 
description of this class of critical spin chains. The ex- 
act AKLT state 25| for a special version of the 5 = 1 
chain was similarly important in confirming Haldanc's 
conjecture of the qualitative differences between half-odd 
integer and integer spins Exact solutions are rare, 
however, and mostly limited to ID models. 

Another, more general approach is to study Hamilto- 
nians numerically. Exact diagonalization of the Hamil- 
tonian is possible only for very small systems; currently 
up to 42 S = 1/2 spins [H [27|. White's density ma- 
trix renormalization (DMRG) scheme [28| and related 
approaches based on matrix-product states [2^] have en- 
abled more detailed studies of the low-energy physics of 
ID systems, including also "ladders" of several coupled 
chains [30]. QMC studies are also competitive in ID in 



3 



the absence of sign problems |3lM3c 

While there has been some progress for 2D systems 
with DMRG 



34[ and tensor-product [35M37| states, 



which arc higher-dimensional generalizations of the 
matrix-product states generated in DMRG [HI, in gen- 
eral these methods have not yet reached the point where 
one can obtain unbiased results for generic nontrivial 
models. There are encouraging developments, however, 
of calculations addressing challenging systems and ap- 
proaching the level of accuracy where definite conclusions 
can be drawn [H,|40[. 

Presently, however, the only numerical approach with 
which one can routinely reach sufficiently large 2D and 
3D lattices in a completely unbiased way is QMC simu- 
lations [4lL [4^| : but, as already noted, they are restricted 
to systems free from sign problems. The class of mod- 
els for which the sign problem is either absent or evad- 
able still contains a vast range of Hamiltonians with non- 
trivial and interesting ground states and quantum phase 
transitions. To study low-energy emergent properties, 
the interactions do not necessarily have to correspond 
in detail to any particular real material — although some 
times that is also possible [43|. The idea is to design a 
sign-problem free Hamiltonian in such a way that it con- 
tains a particular macroscopic (low-energy) phenomenon 
of interest — so that universal physics is captured. 

Bench-mark results obtained from such designer 
Hamiltonians, representing various physical phenomena 
in a prototypical manner, can be very useful experimen- 
tally. Here our main focus will be on the theory side. 
With unbiased large-scale QMC simulations, one can test 
field theories proposed to capture specific classes of quan- 
tum many-body states or quantum phase transitions. 
Moreover, explorations of designer Hamiltonians can also 
serve as "experiments" for discovering novel phenomena, 
and, thereby, stimulate further theoretical developments. 

A designer Hamiltonian corresponds to an effective 
classical statistical-mechanics problem (with a real ac- 
tion) , which arises out of the construction of the sampling 
space in the QMC scheme implemented. In some cases 
quantum phenomena in d dimensions can be studied us- 
ing Monte Carlo simulations with the simplest general- 
ization to a classical model with the same global sym- 
metries in d + 1 dimensions [3 - 14(1 ] (e.g., the thermal 
phase transition of the 2D classical Ising model is in the 
same universality class as the quantum phase transition 
of the transverse-field Ising chain 0,113). Often, how- 
ever, the effective classical model resulting from the con- 
struction of a QMC configuration space, in a complete 
treatment of a given a quantum Hamiltonian, has un- 
usual degrees of freedom and interactions, unanticipated 
in the study of classical statistical mechanics. The naive 
mapping with just d — > d + 1 can then miss important 
and intriguing phenomena special to quantum mechanics 
[22I l23l liq . Thus, in our view the approach of directly 
studying quantum mechanical designer Hamiltonians is 
an important direction in forming an unbiased picture of 
quantum criticality and emergent phenomena. 



B. Outline of the Review 

In this Review, we will focus on some recent exam- 
ples where it has been possible to study quantum phase 
transitions in detail by QMC simulations and to make 
quantitative comparisons of extracted universal quanti- 
ties with field theoretic results. We also discuss other 
fascinating quantum phenomena, e.g., emergent gauge 
fields and entanglement, where QMC studies are playing 
a critical role in establishing benchmarks and gaining in- 
sights that are hard or impossible to extract from the 
field theories. 

In Sec. [H] wc discuss the sign problem and some of 
the QMC schemes used in large-scale studies of quantum 
spin systems, including treatments of S\J(N) symmetry 
and recent progress in computing the Rcnyi entanglement 
entropy. In Sec. IIIII we give examples of conventional 
and deconfined quantum-critical points in SU(2) invari- 
ant models, and in Sec. IIVI we show how generalization 
to SU(iV) symmetry can be used to directly connect to 
large- N studies within field theory. In Scc.|V]wc discuss 
spin liquid states and associated quantum phase transi- 
tions in U(l) symmetric, highly frustrated "XXZ" spin 
models (or, equivalently, interacting hard-core bosons). 
We conclude in Sec. IVII with some further discussion and 
an outlook on future prospects and challenges. 



II. MODERN QMC METHODS 

There is a wide range of QMC methods available for 
studies of different classes of lattice models |U HH, |49[ . 
In this Review we focus on spin models, and therefore 
restrict the discussion to the types of methods most suit- 
able for them. These methods are also often well suited 
for studies of bosons, where a substantial body of work 
has been carried out in the past several years, stimu- 
lated by ultracold atoms in optical lattices [50l - F53| . For 
fermions, there are alternative approaches, based upon 
auxiliary fields decoupling the interactions and subse- 
quent analytical tracing-out of the fermions [4!| |54| . 
These methods have a much less severe fermion sign 
problem (in some cases completely avoiding it, e.g., for 
particle-hole symmetric systems) — with the exception of 
ID systems, where the methods discussed here also work 
well (since the signs from anticommutation normally do 
not appear in this geometry). We will only give a brief 
non-techical overview of QMC methods here, with refer- 
ences to more detailed accounts. 



A. The sign problem 

Solving many-body quantum mechanical problems on 
a computer is generally exponentially hard (in comput- 
ing time versus system size) with known algorithms. The 
statistical mechanics of a d-dimensional quantum Hamil- 
tonian can be rewritten in terms of some classical vari- 
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ables C in d + 1 dimensions, i.e., the partition function 
(T > 0) or the normalization (T = 0) can be cast into 
the form 

z = Y J w c . (1) 

c 

Systems of interest can, for practical, known ways of 
transforming into this form, have both positive and neg- 
ative weights Wc, invalidating the usual interpretation 
of Wc/Z as a probability distribution in Monte Carlo 
simulations. Formally, this difficulty can be side-stepped 
by not including the sign in the probability, using Pc oc 
| Wc | , and compensating for this by weighting observables 
with the sign; 

(O) = £c Wc ° c = Ec°cSign(Wc)\W c \ 
{ ' Ec w c Ec Si g n ( w c) I W c \ 

= (Sign - Q)\ W \ 

(Sign) m ■ [ ' 

The problem here is that the numerator and denomi- 
nator in the last expression both approach zero expo- 
nentially in system size and inverse temperature, and, 
hence, resolving the ratio within statistical noise is expo- 
nentially hard. The advantage of importance sampling 
(time speed-up from exponential to polynomial scaling) 
is then lost. This predicament is the QMC sign problem 

Here, we consider classes of quantum Hamiltonians for 
which it is possible to choose a basis in which all the 
Wc > 0. These sign-problem free models can generally 
be simulated in polynomial time by importance sampling, 
allowing access to large system sizes (with the excep- 
tion of "glassy" systems and other systems where even 
classical Monte Carlo simulations scale exponentially). 
In some cases, sign problems that appear intractable at 
first-sight can be solved in more sophisticated ways, e.g., 
in the Mcron algorithm (55j . Here we implicitly consider 
cases where the "de-signing" is essentially automatic, be- 
ing a direct consequence of the non-positivity of the off- 
diagonal matrix elements of the Hamiltonian. 

B. T > and T = QMC methods 

In finite-temperature methods, the goal is to compute 
thermal averages 

(A) = ~Tr{Ae-e B }, Z = Tr{e^ H } > (3) 

where [3 = 1/T. As an alternative to studying the ground 
state in the limit of T — > 0, in a ground-state projector 
method some operator P(/3) is applied to a "trial state" 
\^o), sucn that \typ) = P(f3)\^ ) approaches the ground 
state when (3 — > oo. An expectation value 

{A) = ±{Hf /} \A\* p ), Z=(*p\*p), (4) 



tends to the true ground state expectation value, (A) —> 
(0|A|0). For the projector, one can use the imaginary- 
time evolution operator P(f3) ~ c~^ H or a high power of 
the Hamiltonian; P(m) = H m . Here m cx (3N gives the 
same rate of convergence for the two choices of projectors 
for a given system size N. This follows from a Taylor 
expansion of the time evolution operator, which for large 
/3 is dominated by powers of the order n = P\Eq\, where 
Eq is the ground state energy. 



C. Path integrals and stochastic series expansions 

The task is now to rewrite Z in Eqs. ([3]) or (U) in 
the form ([I]), and expectation values as Eq. ([2]), with 
classical variables without diagonalizing H. In practice, 
T = and T > schemes for a given model are often 
very similar. In both cases, the exponential operator 
can be treated with path-integral methods — world lines 
in discrete 0, HH or continuous 0, H3, [H[ imaginary 
time — starting from a product of "time slice" evolution 
operators, 

^=fie-- A = A. (5 ) 

Complete sets of states are then inserted between the 
exponentials (and the trace is further taken at T > 0). 
For small A, the matrix elements of the exponentials can 
be evaluated approximately, giving the form ([lj with Wc 
a product of M matrix elements. In modern methods the 
limit A — > is taken at the algorithmic level [1, [ill H3 
and the Monte Carlo sampling is of paths in continuous 
imaginary time. Older approaches employed the Suzuki- 
Trotter decomposition [47| in Eq. ([5]), which typically 
leads to an error C(A 2 ). 

An alternative to Eq. ([5]) is to start from a series ex- 
pansion @, 

°° on 

e-? H = J2?-(-Hr, (6) 

n=0 

and insert complete sets of states between each instance 
of H. The power n itself is importance-sampled along 
with the terms of H (Stochastic Series Expansion; SSE). 
This expansion is sharply peaked around (n) = ^j3E 1 
where E is the total internal energy (under the assump- 
tion that there is no sign problem) , which leads to a com- 
putational effort scaling as @N , which is the same as in 
the path-integral approaches. Again, Wc in Eq. ([XJ) is a 
product of matrix elements. 

The continuous-time path integral can also be seen as 
a variant of SSE, where the evolution operator is writ- 
ten in the interaction picture, with H = Hq + V and 
expanding only in the perturbation V to some conve- 
niently chosen H [HEI. The full SSE, where H Q = 0, 
V = H , is normally more efficient for spin systems, while 
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the continuous-time variant should be better for certain 



boson systems [61 



Regardless of how the exponential operator is treated, 
the difference between Eqs. ^ and (0]) is essentially only 
in the boundary conditions in the imaginary time direc- 
tion; periodic for T > (due to the trace being taken) 
and dictated by the nature of the trial state in T = 
projections. An equal superposition over all basis states 
corresponds to fully open boundaries, while other choices 
of the trial state leads to "biased" boundaries with some 
weighting of the bra and ket states (depending on the 
details of the trial state). The projector approach with 
P{P) = c~ flH is often called the Path Integral Ground 
State method (PIGS) in the context of continuous-space 
systems foU [o3j. 

There is no sign problem in these approaches for spin 
systems with non-positive definite off-diagonal matrix el- 
ements (for a properly chosen trial state in the case of 
the T = approach) . For bipartite antiferromagnetic in- 
teractions a sublattice rotation can accomplish this, or, 
cquivalently, one can note that for such models the signs 
are "invisible", because all non- vanishing terms of the 
path integral or series expansion have an even number of 
off-diagonal matrix elements. For further details of sign 
issues, see, e.g., Refs. [El, [H Ej] ■ Note that there are no 
restrictions on the diagonal interactions in the chosen ba- 
sis, so that one can study, e.g., spin and boson systems 
in which the potential energy is highly frustrated (as we 
will discuss below in Sec. W} . 



D. Sampling with loop and cluster algorithms 

A breakthrough leading to today's efficient QMC al- 
gorithms was the realization by Evertz et al. j8[ that 
the classical Swendsen-Wang (SW) cluster Monte Carlo 
method 64| could be generalized to loop algorithms (i.e., 
generating clusters in the form of loops) for vertex mod- 
els, and that these models are very similar to the effective 
statistical mechanics problems arising in the transforma- 
tions of quantum problems discussed above. While the 
loop algorithm has been adapted to many different mod- 
els^ both within the path-integral [|, El, [§1 ^ and SSE 
jlll | frameworks, they cannot be applied to all cases (in 
analogy with the SW method, which also has limited 
applicability). The loop concept has also been further 
generalized, however, to "worms" [l(| H3, [H[ and "di- 
rected loops" [H],|67|, which are essentially loops that can 
self-intersect during their construction and incorporate 
the detailed-balance principle in a more general manner. 
Recently such a n alg orithm was also formulated without 
detailed balance [68( . In addition, for some models where 
loop algorithms and their generalizations are not appli- 
cable (e.g., transverse- field Ising models), other efficient 
generalizations of the SW cluster concept have been de- 
veloped mill. 



E. SU(2)-symmetric models 

A particularly convenient class of trial states in pro- 
jector QMC studies of SU(2) invariant interactions is the 
amplitude-product states in the overcomplete valence- 
bond (singlet pair) basis (7lj . which were recently gen- 
eralized to include bond correlations [72|. These total- 
singlet states incorporate Marshall's sign rule in a con- 
venient way and also automatically have the appropri- 
ate momentum for the ground state of a periodic system 
(thereby from the outset filtering out a significant frac- 
tion of the excited states). The Hamiltonian in this case 
can typically be expressed in terms of two-spin singlet- 
projector operators (individual ones or products of two 
or more of them) ; 



Pi, 



j)(Si 



(7) 



where Sij denotes the singlet state of spins i and j. A 
power of the Hamiltonian H m , for fixed m or in a series 
expansion, Eq. © , is expanded into all possible products 
of these projectors, which act on the ket trial state in (QJ. 
In early versions of this scheme, the resul ting pa ths of va 
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lence bond states were sampled directly 731475 1. while in 
more recent formulations the diagonal and off-diagonal 
parts of the projectors are also sampled, along with spin 
configurations compatible with the valence-bond configu- 
rations of the trial state (which are also sampled) [76| . An 
example of a configuration in such a simulation is shown 
in Fig. [TJ This configuration space is amenable to the 
very efficient loop updates discussed above [l2|. Apart 
from the boundary conditions in the "propagation" direc- 
tion (and the associated sampling of the valence bonds 
in the trial state), this kind of projector method is in 
practice implemented in a very similar way as the SSE 
method for SU(2) models (42J . 

Evaluation of expectation values ( "measurements" ) in 
projector QMC simulations are carried out at the mid- 
point of the configuration (while in T > methods aver- 
ages can be taken over the whole time range, due to the 
time-periodicity), as indicated in Fig. [TJ In the valence- 
bond basis, when propagating the bra and ket states in 
the pure valence-bond basis (which corresponds to av- 
eraging over all compatible spin configurations), most 
quantities of interest can be expressed using the loops 
of the transition graph formed when superimposing the 
projected bra and ket bond configurations [7l|, [77[ . Es- 
timators for operators A expressed in the space {£} of 
transition-graph loops are typically of the form 



A(C) 



(Vl\A\V r ) 

(Vi\V r ) 



(8) 



where the overlap of the two valence-bond states is given 
by 



(V l \V r )=2 N °- N '\ 



(9) 



where iV D is the number of loops and N the system size. 
The numerator in Eq. JSJ depends on the loop structure, 
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FIG. 1. A configuration in a valence-bond projector simula- 
tion of a 4-site Heisenberg chain; here with a small projec- 
tion power, m = 2 [l^]. The vertical bars indicate singlet- 
projectors operators, while the arcs on the left (bra) and 
right (ket) represent valence bond configurations of the sam- 
pled trial state. Open and closed circles indicate up and 
down spins, which are unchanged between operators and rep- 
resented by horizontal lines there. These lines, along with 
the valence bonds, form loops that can be nipped (i.e., all 
spins traversed by the loop are flipped) without changing the 
configuration weight. A configuration in a T > simulation 
corresponds to imposing periodic boundary conditions in the 
"time" direction (reflecting the trace operation), instead of 
the independent boundaries terminated by valence bonds. 



e.g., for A — Si ■ Sj, the estimator A(C) is if sites i 
and j belong to different loops and ±3/4 for sites in the 
same loop (with + and — for the same and different sub- 
lattices, respectively). Hig her-order correlation functions 
are discussed in Refs. |78| and (79j . 



F. Generalizations from SU(2) to SU(N) 

The above methods for SU(2) models can be extended 
to a certain class of spin models with an SU(A^) sym- 
metry on bipartite lattices. First, note again that the 
SU(2) Heisenberg operator is simply a singlet projector, 
Eq. ((7j). In order to generalize this operator to SU(iV), 
we choose models that have an SU(iV) spin with a funda- 
mental representation on the A sub-lattice and an SU(iV) 
spin with a conjugate to the fundamental representation 
on the B sub- lattice [8(J HJ • Denoting the N states on 
each site by \a), where 1 < a < N, and the generators of 
the fundamental representation by T a (N x N matrices), 
the SU(7V) singlet of a spin i on the A sub-lattice and j 
on the B sub-lattice is given by, 



5Zk 



(10) 



We note here that for N = 2 and on bipartite lattices, 
this singlet is equivalent to the standard way of writing 
the singlet as (| fj,} — | J,f))/v2, by making the transfor- 
mation | — >• | I), | I) —> — | t) on onc °f the sublattices. 
We can use \sfj) to rewrite the manifestly SU(iV) invari- 
ant version of Eq. (J7J) as an explicitly sign-problem free 
projection operator, 



^ T a. T *. 



1 

]V2 



p, 



y ■ 



(11) 



generalizing the Heisenberg interaction for two spins on 
opposite sublattices. 

Another interaction of interest is an SU(iV) invariant 
interaction between sites having the same representation; 
the permutation operator, which is defined in the follow- 
ing way by its action on a ket: 



—UiAotiPj) = \PiOtj). 



(12) 



We can now relate this operator to the generators T a : 

1 



\ rpa rpa 

2-^1 » ' 3 



N 2 ~ Uij - 



(13) 



This is an SU(A r ) invariant generalization of the SU(2) 
Heisenberg interaction for two spins on the same sublat- 
tice. 

Hamiltonians based the SU(7V) operators and Hy 
are sign-problem free if they come with negative signs, in 
which case they generalize the SU(2) antiferromagnctic 
and ferromagnetic Heisenberg exchange, respectively. In 
terms of the definitions of Py and n^, the algorithms 
presented in Sec. Ill Bl can be simply generalized from 
SU(2) to SU(AT) by extending the number of "colors" 
(of the spins and the loops) from 2 to iV (82W8i |. As 
a consequence of this, estimators expressed in terms of 
transition-graph loops are also modified, e.g., in the state 
overlap Eq. © 2 is replaced by N. One can even gener- 
alize simulations in the pure valence-bond basis to non- 
integer N H. 



G. Renyi entropies via the replica trick 

In addition to conventional physical observables eval- 
uated according to Eqs. (J3]) and (g]), e.g., various corre- 
lation functions, QMC methods have recently been de- 
veloped that are capable of measu ring the degree of en- 
tanglement in a quantum system [851 ] - Several quanti- 
ties related to entanglement have been explored in the 
QMC context for their abilities to identify and charac- 
terize quantum phases and phase transitions [86l - [9lT |. In 
this Review, we concentrate of measures of the entangle- 
ment entropy, specifically the Renyi entropies [92| , 



I- a 



hi 



(14) 



for integer a > 1. Here, pA is the reduced density matrix, 
Pa = Tr b{p}, and A is a subregion of a lattice system 
(with B being its complement). The von Neumann en- 
tanglement entropy corresponds to the limit a — > 1. 

The direct measure of Renyi entropies can not be done 
in QMC using conventional estimators, however recent 
work has demonstrated that it is possible to measure S a 
for a > 2 using a replica trick [93j-|97[. Here, a copies 
of the simulation cell arc used, with modified periodic 
boundaries in the "propagation" (or imaginary time) di- 
rection. The evaluation of Renyi entropies proceeds in 
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formally different ways for T = and T > QMC meth- 
ods hR [la ]. Namely, for T = 0, the calculation of S a 
proceeds in analogy with Eq. |4]) (where is the wave- 
function of the a-times replicated system), with the op- 
erator A replaced by a "swap" operator for a = 2 [lj| or 
permutation operator for a > 3 [981 ] - 

The actual simulation procedures are easiest to under- 
stand in the context of the projector algorithm discussed 
in Sec. Ill El Here, the second Renyi entropy (for exam- 
ple) is given by the mid-point evaluation of a SWAP^ 
operator on two replicas of the system (see Fig. that 
literally swaps states (spins or valence-bond endpoints) 
between replicas when they lie in region A. Then, 



5 2 = -ln (SWAP^(£)) , 



(15) 



where the expectation value is calculated at the mid- 
point of the configuration, and expressed using the loops 
of the transition graph C formed when superimposing the 
projected bra and ket states; Eq. ©. Namely, 



(Vi\SWAP A \V r ) 
(Vi\V r ) 



O-Wswap— N 



(16) 



where A^ swap is the number of transition-graph loops of 
the swapped configuration, while A is for the configura- 
tion before the swap operator is applied (Fig. [2J). We re- 
fer to Ref. [98| for further details of this particular T = 
projector QMC implementation. 

At finite temperature, the explicit evaluation of the 
SWAP^ operator can be replaced by an evaluation of 
the difference in free energies between a replicated sys- 
tem (similar to Fig. [2]), and an un-replicated system [lq |. 
using Eq. ([3]). Although this technique looks quite dif- 
ferent formally, both the T — and T > techniques 
can be viewed on a similar footing when considered in 
the context of SW loop or cluster algorithms, much like 
the unified framework of the QMC itself, discussed above 
in Sec. Ill Bl In some cases, the implementation for these 
two different techniques can be made to be practically 
identical. 



III. GROUND STATES AND QUANTUM 
PHASE TRANSITIONS IN SU(2) MODELS 

It has been a matter of debate for a long time whether 
the non-linear cr-model can correctly capture a phase 
transition from the Neel state to a non-magnetic state 
in two dimensions for 5 = 1/2 spins 20]. Other terms, 
related to Berry phases of the spins, have to be included 
in t he lo ng wavelength description to obtain a VBS state 
[99L llOOj . It was recently proposed that the Neel and VBS 
states are separated by a "deconfined" quantum-critical 
(DQC) poin t [22"! described by the non-compact CP 1 field 
theory \l0l\ (in contrast to the nonlinear cr-model, which 
is equivalent to the compact CP 1 description, i.e., in the 
deconfined scenario the Berry phases render the gauge 
field effectively non-compact at the critical point). This 
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FIG. 2. The SWAPa operator for calculating the Renyi en- 
tropy S2, where subsystem A consists of the two top sites. The 
simulated system is composed of two non-interacting replicas, 
top (compare to Fig. [T]), and bottom, which are sampled in- 
dependently. Just to the right of the mid-point of the con- 
figuration (dashed line) the SWAPa operator (green vertical 
bar) connects the A spins in the top replica to the bottom 
one, as indicated by the colors of the loops [9^ . 



field theory describes two flavors of spinons interacting 
with a non-compact U(l) gauge field. The VBS and Neel 
order parameter form due to confinement and condensa- 
tion, respectively, of spinons, and therefore cannot be 
regarded as two separate order parameters at the phase 
transition. 

In this section we first discuss the rather well-studied 
continuous quantum phase transition between the Neel 
state and a quantum paramagnet in systems with bi- 
modal Hcisenberg couplings forming a static pattern. In 
this case the paramagnet does not break any symmetries 
of the Hamiltonian and the transition is believed to fall in 
the 3D 0(3) universality class. A VBS state with sponta- 
neously broken lattice symmetries can be achieved with 
frustrated interactions, which, however, are si gn p rob- 
lcmatic in QMC simulations. The J-Q model [102| is a 
designer Hamiltonian circumventing this problem. We 
will review some of the recent work aimed at character- 
izing its VBS state and Neel- VBS transition. 



A. Neel— paramagnetic transition in dimerized 
systems 

Consider 5=1/2 spins on a 2D bipartite lattice (e.g., 
the simple square lattice) , interacting with antiferromag- 
netic Heisenberg exchange of strength J\ , written using 
the singlet projectors (JT]), 



H 



- j i Pi * 



(17) 
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Here denotes nearest-neighbor spins. This system 

has long-range Neel order at T = 0. Now introduce 
dimerization, by considering pairs (dimers) of spins 
such that each spin belongs exactly to one dimer. In ad- 
dition to the nearest- neighbor couplings (fT"7|). let there 
be additional intra-dimer couplings of strength A (i.e., 
the total intra-dimer coupling is J\ + A), such that the 
Hamiltonian is 



H 



(18) 



If A > and J\ — > 0, the ground state is clearly a prod- 
uct of singlets on the dimers — a quantum paramagnetic 
state that we will refer to as a valence-bond-liquid (VBL), 
motivated by the fact that the state breaks no symme- 
tries and is strongly interacting away from the large-A 
limit. When A = the system is the standard Hciscnberg 
model (|17p . There is, thus, a quantum phase transition 
between Neel and VBL ground states as a function of 
g = (Ji + A)/Ji. The dimers can be arranged in differ- 
ent ways, and instead of dimers one can also use larger 
Ji-coupled units of an even number of spins, e.g., 4-site 
plaqucttes, or one can use a bilayer, with J\ and Ji the 
inter-and intra-layer couplings, respectively. 



ing behavior in this extended T > region is an impor- 
tant generic characteristic of quantum phase transitions, 
in contrast to classical phase transitions with typically 
very narrow regions of criticality. This aspect of dimer- 
ized mo dels has also been investigated in detail in QMC 
studies dTMlJ 



3. Experimental realizations 

A Neel- VBL transition of the type discussed above can 
in principle be experimentally realized in systems con- 
sisting of coupled dimers as a function of pressure. The 
authors are not aware of any material where the phase 
transition can be crossed in a quasi-2D systems of weakly 
coup l ed p lanes, but a well-studied 3D case is TICUCI3 
[l!5l Ill6l |. Here neutron scattering shows a transition 
from a paramagnet to a Neel state at a critical pressure 
[l!7j . including the analogue of the Higgs boson (lon- 
gitudinal mode) in the ordered state [l!8| . In this case 
the values of the couplings and their pressure dependence 
are not known in detail, but one can still learn a lot from 
studies of designer H amiltonians by focusing on universal 
properties 



^sign er f l; 



1. T = criticality 

Based on symmetry arguments alone, the Neel- VBL 
transition in these dimerized (or polymerized) systems 
should be in the classical 3D 0(3) universality class. 
Many QMC studies have been devoted to checking the 
expone nts again st available results for the classical tran- 
sition 103l4l08| . For some of the systems, e.g., bi- 
layers 106l and columnar dimers on the square lattice 
[lallOa 1081 ] , results have been obtained that rival classi- 



cal calculations in precision (statistical error bars) and in 
general the exponents agree very well with each other. In 
other cases, e.g., with the dimers arranged in a staggered 
pattern, som e studies initially indicated a new universal- 
ity c lass |l07| . while others supported the 0(3) transition 
[I09j . Later it was realized that these cases, where the 
dimer pattern lacks inversion symmetry, are associated 
with certain operators (cubic interactions) in the field- 
theory formulation which are not present in systems with 
inversion symmetry [llCf . These operators are asymptot- 
ically irrelevant in the RG sense, but can lead to substan- 
tial finite-size corrections that can easily be mistaken for 
a different universality class. 



B. Deconfined Neel— VBS transition in J-Q models 

While spontaneously formed VBS states in transla- 
tionally invariant 2D spin syste ms have b een discussed 
for more than two decades [8l|, 1 1241 . Il25|, most of the 
early work focused on frustrated systems 125l4l27 |. for 
which it is very difficult to carry out unbiased numerical 
calculations. The proposal of a DQC point separating 
the Neel and VBS states motivated detailed numerical 
studies of this transition, and, thus, construction of suit- 
able designer Hamiltonians amenable to QM C sim ula- 
tions. Numerical work preceding the theory [l28j | had 
already indicated the intriguing possibility of continuous 
magnetic- VBS quantum phase transitions (where stan- 
dard arguments based on Landau- Ginzburg theory pre- 
dict generically first-order transitions) in U(l) (quantum 
XY) models. Tractable SU(2) spin models exhibiting 
Nccl-VBS transitions were lacking, how ever, until the 
introduction of the J-Q class of models [102], where J 
refers to Hciscnberg exchange on a bipartite lattice and 
Q is a multi-spin interaction which competes against Neel 
order but does not generate a sign problem. 

The Q interactions are most conveniently expressed in 
terms of the singlet projectors, Eq. ([7]), in the simplest 
case on the 2D square lattice; 



2. T > scaling 

Field-theoret ic d escriptions using the nonlinear a- 
model [l9|, lllll ] have given a wealth of predictions 
also for the T > quantum-critical "fan" extending out 
into the plane (g,T) from the critical point (<7 C ,0). Seal- 



(ijkl) 



(19) 



where the site pairs ij and kl form parallel edges (with 
both horizontal and vertical arrangements included) on 
a 2 x 2 plaquctte. One can also consider more than two 
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FIG. 3. Q2 and Qs interactions on the square lattice. The 
bars indicate the locations of singlet projectors dj on site 
pairs ij. The Q2 and Qs arrangements in the left and center 
panels lead to columnar VBS states, while the right Qz ar- 
rangement leads to a staggered VBS. The Hamiltonian con- 
tains all unique translations and 90° rotations (and in the 
staggered case also reflections) of the operator patterns. 



projectors in the product and denote by Q p an interaction 
with p projectors. Fig. [3J illustrates Q2 and Q3 terms 
leading to columnar VBS ground states. In addition to 
the number of projectors in the product, they way they 
are arranged relative to each other on the lattice is also 
crucial, with different patterns leading to different VBS 
ground states of the "pure Q" models (fl~9|) . VBS to Neel 
transitions can be studied in J-Q models, H = Hj + Hq, 
as a function of the ratio J/Q of the standard Heisenberg 
exchange and the multi-spin interaction. 

A staggered Q3 term (illustrated in F ig. [3]) on the 
square lattice was investigated in Ref. [l29l | . This leads to 
a strongly first-order transition. In contrast, the colum- 
nar Q 2 and Q3 arrangements in Fig. [3] lead to continuous 
transitions (or, in principle, the transition could be very 
weakly first-order). Various types of Q term s hav e also 
been investigated in the honeycomb lattice [l30| . with 
sim ilar r esul ts. In analogy with quantum dimer mod- 
els jl3ll . 1 13 j |. first-order transitions correspond to VBS 
states that do not support any local dimer fluctuations, 
while the continuous transitions are into VBSs where 
such fluctuations exist. Alternatively, one can relate the 
different behavio rs to diffe rent types of topological de- 
fects in the VBS pol . fl33] . 



1. T = Critical behavior 

The critical behavior of both the J-Q2 and J-Q3 mod- 
els (with the columnar dimer arrangements in Fig. [3]) on 
the simple 2D square lattice has been analy zed in detail 
in several different ways [zl, [To3, HMCM- 0ne useful 
quantity for studying the destruction of the Neel state is 
the spin stiffness p s . It should scale at a critical point 
as p s ~ L 2 ~ d ~ z , where z is the dynamic exponent. The 
DQC theory is Lorenz-invariant, i.e., z = 1. One can test 
this prediction by graphing Lp s (L) versus the coupling 
ratio J/Q for different Lx L lattices (in the ground state 
or with the inverse temperature /3 oc L z ). Such curves 
should cross at the critical point. In pr actic e, crossing 
points often exhibit some drift with L [l36| . reflecting 
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FIG. 4. Finite-size scaling of the spin stiffness of the J- 
Q2 model in the vicinity of its quantum-critical point; from 
Ref. [7^|. The calculations were done at inverse temperature 
/3 — L. The standard critical scaling form Lp s ~ constant 
has been modified by a log-correction (with Lo = 0.9) in or- 
der to compensate for a weak drift of the crossing points. 
The red curves, which cross at a single point, the critical 
point q c = (J/Q) c = 0.0447(2), represent a common scaling 
function for L > 48 (approximated by a common polynomial 
fitted to all the data points), f[(q — q c )L x l v \, with v = 0.59. 



scaling corrections [76|]. Results for the J-Q2 model are 
graphed in Fig. @] and discussed next. 

In the J-Q models the Lp s -cro ssing drift is anoma- 
lously large, which in one study [134{ had been inter- 
preted as a weakly first-order transition (for which Lp s 
should eventually diverge as L, due to coexistence of a 
stiff and a non-stiff phase). Another possibility is that 
the model has unusually large scaling corrections, that 
could be intimately associated with the nature of the 
DQC point. For example, in the DQC theory the per- 
turbation to the critical point leading to VBS order (a 
quadrupled monopole operator) is dangerously irrelevant. 
It could potentially lead to large, but ultimately conven- 
tional scaling corrections. Another possibil ity is th at the 
corrections are logarithmic (multiplicative) [7a , ll38t | , sim - 
ilar to those known in the critical Heisenberg chain [l39j | . 
Such corrections do not appear in the large-iV calcula- 
tions of the CP^ -1 theory, but potentially they could 
appear for small N, as has recently also been argu ed for 
based on a modified version of the DQC theory [140| . Nu- 
merical data for the J-Q2 model [7fj, on L x L lattices 
with L up to 256, can be well accounted for by a log- 
correction, as shown in Fig. 21 but the deviations from 
pure scaling can also be fitted with a conventional mul- 
tiplicative correction of the form (1 + aL^^) with large 
a and small lo. 



2. T > critical scaling 

The quantum-critical "fan" in the (J/Q,T) plane has 
also been investigated. Initial calculations confirmed that 
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the scaling is govern ed by a z — 1 critical point [1361 ]. 
Later studies [76l . lll4 [ noted log-like corrections also here, 
e.g., while the conventional z = 1 scaling implies that 
the uniform magnetic susceptibility scales linearly in T 
[2p| , the behavior actually appears to be of the form \ ~ 
T[l + aln(l/T)]. These corrections, like those at T = 
discussed above, remain an intriguing unexplained aspect 
of the Neel-VBS transition in the J-Q models. 



3. Emergent U(l) gauge fluctuations 

An important and remarkable aspect of the DQC the- 
ory is the emergent conservation law of the gauge flux, 
due to the irrelevance of monopolcs at the critical point . 
There is a U(l) symmetry associated with this conserva- 
tion law at the critical point, which is also emergent and 
not explicitly contained in the J-Q model, which only has 
an SU(2) spin-rotational symmetry, the discrete symme- 
tries of the square lattice and time reversal. The emer- 
gent U(l) symmetry can, however, be o bserv ed in the 
fluctuations of the VBS order parameter [l02j |. Naively 
one would expect these fluctuations to reflect the Z4 sym- 
metry of the VBS (which is dictated by the square lat- 
tice), but instead the symmetry group expands to U(l) 
upon approach to the critical point. The interpretation 
of this emergent U(l) VBS symmetry is that it reflects 
the non-compactness of the gauge field ( "photon" ) of the 
DQC theory. Here we will review how this emergence is 
detected in QMC simulations. 

Consider the VBS (dimer-dimer) correlation function, 



C a (r tJ ) = (B a (Ti)B a (Tj)), 



(20) 



where = Yi — Vj is the spatial separation of the oper- 
ators B ai a = x,y, which measure the spin correlations 
on nearest-neighbor bonds in the a direction; 

£*(r) = S(r)-S(r + x), B fi (r) = S(r) -S(r + y). (21) 

We define the expectation value of this operator when 
averaged over all sites r = (x, y); b a = (B a (x,y)) 7^ 0. 
In a VBS, the expectation value acquires a spatial mod- 
ulation. With the dimers (stronger bonds) oriented in 
the x direction, the asymptotic form of the spatially 
averaged x dimer correlations in a columnar VBS is 
C x {x,y) = bl + {-l) x {±D 2 +ae- r ^), where r = (x 2 +y 2 ), 
D is the magnitude of the order parameter, and a is a 
constant. The rate of decay of the oscillating term to- 
ward a non-zero value as r — > 00 defines the correlation 
length £ (where we suppress the directional dependence 
of £, which can be parameterized by £ x and £ y ). This 
is the standard definition of a correlation length, as the 
length-scale associated with the order parameter. Nor- 
mally this (the larger of £ x ,£y) is the largest length-scale 
in the system. 

In a symmetry-broken x-oriented VBS, we can consider 
also dimer correlations in the direction y perpendicular to 
the dimers. The asymptotic form of these correlations is 



Cy(x, y) = 6| + fce _r / A , where k is another constant (and 
again there are really two decay constants; A x and A y ). 
These correlations are, at first sight, not associated with 
the order parameter, and the length-scale A would then 
normally not be referred to as "the correlation length" . 
However, A can actually be larger than £, and this is 
at the heart of the emergent U(l) symmetry of the VBS. 
According to the DQC theory, both lengths diverge as the 
critical point is approached, and A ~ £ 1 + a with a > 0. 
Close to the critical point, where the order parameter 
D is small and A 3> £, on length scales less than A the 
system will appear to have both x and y VBS order. 

One can also associate A with fluctuations of the angle 
of the order parameter. Since the VBS order parameter 
is a vector, D = (D x ,D y ), one can define a magnitude 
D and an angle 0. A defines a length scale below which 
the order parameter exhibits U(l) symmetry — the sys- 
tem not only has both x and y order up to this length 
scale, but one can think of the angle as fluctuating uni- 
formly, even though it has to take one of the values nn/2 
when the order parameters D x and D y of a symmetry- 
broken state are averaged over regions of size 3> A. There 
are many interesting consequences of the emergent U(l) 
symmetry of t he V BS order parameter, as further dis- 
cussed in Ref. [l4ll | . 

This phenomenon of an emergent symmetry higher 
than that of the order parameter in an ordered state is 
in fact a general aspect of systems with dangerously ir- 
relevant perturbations. These perturbations reduce the 
symmetry of a system in such a way that the standard 
critical exponents at a phase transition into an ordered 
state are unaffected, b ut the sy mmetry broken in the or- 
dered state is reduced |l42l . Il43| . A prototypical example 
is the 3D classical XY model with an added on-site po- 
tential hcos(q8i) on al l sites i, where 0i is the spin angle 
and q an integer [l44j . This model orders at a critical 
temperature T c (h), but, regardless of the value of h, the 
critical exponents are fixed at those of the standard XY 
transition. However, instead of breaking the continuous 
U(l) symmetry, as the original h = model docs, now 
only a q-iold (Z q ) symmetry is broken. As in the VBS 
discussed above, there is a length-scale A ~ £ 1 + a (q) in the 
ordered state, which is associated with U(l) symmetric 
fluctuations of a course-grained order parameter (magne- 
tization) m = (m x ,m y). Finite-size scaling procedures 
have been developed jl45l | to access the exponent a(q) 
directly based on the cross-over from U(l) to Z q symme- 
try seen in the order-parameter distribution P(rn Xl m y ). 
Such techniques have also been employed to study the 
VBS fluctuations of the J-Q model, as we discuss next. 

In simulations of a VBS on a finite lattice, where the Z4 
symmetry is not broken, one can study the emergent U(l) 
and the length-scale A by accumulating the probability 
distribution P(D X , D y ), where D x and D y are defined 
on the whole system in terms of matrix elements of the 
operators B a in Eq. (|21[) . evaluated for individual con- 
figurations using transition-graph loops, as explained in 
Sec. IIII1 and Fourier transforming at (7r,0) and (0,7r) for 
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IV. NEEL-VBS TRANSITION IN SU(AT) SPIN 
MODELS 



FIG. 5. Probability distribution P(D x ,D y ) of the VBS or- 
der parameter of the J-Q3 model on a 32 x 32 lattice at 
J/Q3 = 0.575 in (a) and 0.1765 in (b). The critical ratio 
is (J/Qs) c ~ 0.67. Brighter colo rs correspond to higher prob- 
ability. Adapted from Ref. [l37t ] ■ 



In this section we will discuss quantum criticality in 
SU(7V) spin models on 2D square lattices, as system- 
atic generalizations of the SU(2) systems discussed above. 
Sees. HVA1 ITVBl ITVCl and ITTdI descri be sev eral studies 
carried out on a single layer, and Sec. IIVEI summarizes 
a recent study on a bilayer. 



A. SU(iV) Heisenberg models 



a = x and a = y correlations, respectively. Fig. [5] shows 
resul ts for the J-Q3 model at two different ratios J/Q3 
[l37j | . Far away from the critical point the histogram ex- 
hibits a four-fold symmetry consistent with a columnar 
VBS, while closer to the transition point the distribution 
is ring-shaped. The radius of the ring corresponds to the 
magnitude D of the order parameter, and the absence 
of four-fold symmetry implies that here the system size 
L = 32 < A. By analyzing the U(l)-Z4 cross-over as 
L increases, using finite-size scaling techniques, one can 
extr act t he exponent governing the length-scale A. In 
Ref. [l37| this resulted in A - £ 1+a , with a w 0.2 



4. Spinon deconfinement 

The emergent U(l) length scale discussed above should 
also be reflected in the confinement of spinons (e.g., 
the size of a bound state) in the VBS state. By uti- 
lizing properties of the valence-bond basis, discussed in 
Sec. HIE] and g eneralizations to non-zero total magneti- 
zation jl46l . ll47| . it is possible to directly access 5" = 1 ex- 
citations of VBS states and study their structure in terms 
of bound or unbound (deconfined) spinons. Such calcu- 
lations have already been used to characterize deconfined 
spinons using a de signe r Hamiltonian with a VBS state 
in one dimension [1481 ] . Similar calculations for 2D J- 
Q models should shed further light on the nature of the 
spinons in the vicinity of the DQC point. 

Another approach to confirming spinon excitations, 
which has already been employed for the 2D J-Q model, 
is to study the thermodynamic properties at criticality 
and compare them with predictions for a si mple model 
of a gas of thermally excited bosonic spinons |l!4| . The 
results show a remarkable consistency among different 
quantities which are related through expressions for the 
spinon gas, thus lending strong support to the deconfine- 
ment scenario. 



A natural extension of the SU(2) Heisenberg model is 
to spin models with SU(iV) symmetry. Such generaliza- 
tions were initially introduced to facilitate access to a 
solvable large- A?" limit [10, Hl[. In the meantime, new 
physical systems have also motivated an interes t in such 
models at arbitrary finite values of N [l49l ll5C| , leading 
to extensive studies of SU(A r ) spin systems. The physical 
degrees of freedom of a particular spin model depend not 
only on the global symmetry but also on the chosen rep- 
resentation under which the spins transform. Thus, there 
are a number of distinct extensions of the SU(2) model 
to SU(A), depending on what properties are intended to 
be retained in the generalization. 

In the context of the study of the DQC point at the 
Neel-VBS transition, the ability of two spins to form a 
singlet is the most central characteristic required to form 
a VBS. It is possible to choose SU(A r ) models which have 
this property by working on bipartite lattices and as- 
signing SU(A^) spins that transform as relative conjugate 
representations on each of the sublattices, as discussed 
in Sec. Ill Fl In such models two spins on opposite sublat- 
tices can form a singlet, allowing for the possibility of a 
valence-bond solid. The simplest SU(A^) invariant two- 
site sign-problem free Hamiltonian interaction one can 
define for such models is given by Hi in Eq. (|17|) . with 
Pij generalized to the SU(AT) singlet projector defined in 
Sec. Ill Fl We note again that for N = 2 this reduces to the 
familiar S — 1/2 antifcrromagnetic Heisenberg model. 

The Hamiltonian Eq. (|T7| has been studied extensively 
as a function of N. On a ID chain the well known crit- 
ical Bethe state for N = 2 gi ves way to a VBS ordered 
state for all N > 3 [H Em [H3. The 2D version of 
Hi was shown in the large- A^ limit to map onto a quan- 
tum dimer model with only a kinetic term [l53| |. From 
direct numerical studies of the quantum dimer model, it 
is known that the kinetic-only model orders into a VBS 
[154 Il55l ] . This implies that in Hi there must be a tran- 
sition between the magnetic state at AT = 2 and the VBS 
at N — 00. The location of the transition was first de- 
termined to lie betw een N = 4 and N = 5 by QMC 
simulations (82l . 1 1561 ] . 
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B. SU(JV) Ji-Q model for TV < 4 



In the previous section we saw that on the square lat- 
tice the Neel-order found for TV = 2,3,4 gives way to VBS 
order for TV > 5 in the SU(TV) antiferromagnet, Eq (|T7|) . 
In order to access the quantum phase transition by QMC 
simulations at a given value of TV < 4 one needs new sign- 
problem free "designer" couplings that destroy the Neel 
phase and result in the VBS. The first such term dis- 
covered was the Q-interaction, Eq. (fl~9|) . introduced for 
SU(2) in Sec. IIII Bl It can be generalized for any TV, with 
Pij again the SU(iV) singlet projectors. Since the matrix 
elements of are all positive, such a Q coupling is free 
of the QMC sign problem. 

It turns out that the Q-only model, with the projectors 
arranged as in Fig. [31 is always VBS ordered, so that for 
the cases o f TV = 2, 3, 4, the J\-Q model has a Neel-VBS 
transition [lO^ . [l37l ] . Detailed QMC studies for these TV 
find strong evidence for a continuo us tr a nsition between 
the Neel and VBS ground states [HI EM EH KEf% . 
As discussed above for the SU(2) case, certain observ- 
ablcs show deviations from standard scaling laws for the 
system sizes studied, that could either be due to conven- 
tional corrections to scaling or actual violations of the 
scaling laws (but they do not resemble violations of the 
kind that would be expected at a first-order transition 
[7^1 )■ The origin of these corrections remains to be re- 
solved. 




FIG. 6. Phase Diagram of the SU(TV) symmetric sign- 
problem free J1-J2-Q model as a function of N. Each of the 
couplings has been introduced in the text: the generalized 
nearest-neighbor antiferromagnetic exchange J\ in Sec. II V K\ 
the four-spin coupling Q = Q2 in Sec. HVBl and Fig. H and the 
generalized ferromagnetic next-nearest-neighbor coupling J2 
in Sec. II V CI This SU(TV) antiferromagnet allows for an unbi- 
ased study of deconfined quantum criticality at the Neel-VBS 
transition for each value of TV. 



D. J1-J2-Q model; anomalous dimensions 



C. SU(iV) J1-J2 model for TV > 5 

For TV > 5, the Q term strengthens the VBS that 
is already present in the J\ only model, and hence to 
study the Neel-VBS transition for large-TV, a new de- 
signer Hamiltonian is needed. To supply such a model, 
an interaction J2 can be included between sites on the 



same sublatticc which stabilizes the Neel state 84 



Hi 



E 



n, 



(22) 



where the Hy interaction was introduced in Eqs. (fT2"]) 
and (fl"3|) . In order to see that this interaction favors the 
Neel state, we note that in the H2 model by itself, the 
two sublattices are decoupled. The ground state corre- 
sponds to having an independent SU(TV) ferromagnet on 
each sublattice. Turning on a small J\ <C Ji interac- 
tion will clearly cause the sublattice degeneracy to be 
lifted, resulting in the two independent ferromagnets to 
lock into a single antiferromagnetic state. These argu- 
ments are true independent of the value of TV. Thus, for 
every TV > 5 the VBS state is realized when J\ 3> J2 
and the Neel state is obtained when J 2 3> J\. Thus, 
there must be at least one quantum phase transition as 
the ratio J2/J1 is tuned. From numerical simulations of 
the J1-J2 model, there is compelling evidence for a direct 
transition between the two phases, with no indication of 
an intervening phase 84 j . 



One can also combine the J1-J2 and J-Q model. The 
phase diagram of the resulting J1-J2-Q model as a func- 
tions of TV is summarized in Fig. [5] So far, the phase 
boundaries have only been computed in the two perpen- 
dicular planes shown in the figures. We here discuss some 
aspects of the critical behavior for TV E {2,3,..., 12} 
obtained in these calculations and relate the results to 
large- TV expansions within the DQC theory. 

An important prediction of the DQC scenario [22[ at 
the Neel-VBS transition is that both order parameters 
are simultaneously quantum critical at the phase tran- 
sition and that space and time scale in the same way, 
i.e., the dynamic critical exponent z = 1. This im- 
plies that both correlation functions decay as Lorentz- 
invariant power laws, 



Cjv,v(r,r) 



1 



( r 2 +c 2 r 2)(l+>) V -,jv)/2- 



(23) 



where CV and Cy are the two point correlation func- 
tions of the Neel and VBS order parameters. The indices 
tjn and r\y arc the so-called anomalous dimensions of 
the Neel and VBS order parameters. The DQC scenario 
predicts that the continuum field-theoretic universality 
of the Neel-VBS critical point is described by the non- 
compact CP^ -1 universality. Quantitative estimates for 
the universal indices that characterize this unive rsality 
class arc available only in the large- TV limit 159]. The 
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FIG. 7. Anomalous dimensions of the Neel (left) and VBS 
(right) order parameters as a function of N. The main panels 
show r]M and rjv versus 1/N . For N — 2,3 and 4, the data 
are for the J-Q model [l37l ]. and the results for N > 4 are 
for the J1-J2 model [84l |. The analytic results from the 1/N 
expansion of the CP^ 1 field theory, Eq. (SP, are shown as 
thick red lines. The left and right insets show JV(1 — t/n) and 
(1 + r)v)/N, respectively. These quantities must be finite in 
the N — ¥ 00 limit according to the DQC theory and should be 
given by the prefactors of the iV-dependent terms in Eq. (|24|) . 
The values are indicated by the solid horizontal lines. The 
next corrections to the exponents have not been computed 
analytically yet, but we can estimate them approximately as 
VN = 1 + 32 / (iv 2 N)- 3.6(5) /N 2 ,1 + rjv = 0.24927V + 0.68(4) 
(shown as dashed fitted lines in the insets). Graphs repro- 
duced from Ref. 0. 
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FIG. 8. Physics of the SU(JV) bilayer model, (a) Lattice 
and couplings. White sites form the A sublattice with SU(JV) 
spins that transform in the fundamental representation. Black 
sites on the B sublattice transform as the conjugate to this 
representation, (b) Phase diagram of the SU(6) bilayer model 
as a function of g± — J±/ Ji and g-z = J2/J1. The black circle 
at g± — is the Neel- VBS quantum critical point. Remark- 
ably, in the bilayer geometry with g± 7^ this continuous 
transition becomes first-order (indicated by open circles), due 
to the cancellation of Berry phases between the layers. The 
VBL phase at large g± is a featureless paramagnet with sin- 
glets predominantly forming between the planes. Cartoons of 
the basic singlet physics in the symmetry-broken VBS (left) 
and uniform VBL (right) phases are shown. T he S U(JV) Neel- 
VBL transition is discussed at length in Ref. [l64l [. 



1/N expansions for 77 ^ jl6Cll | and r\y |l00l . Il6l| are, 
32 

^ = 1 -^7> l + r\ v = 2S x N (fcwO.24). (24) 
it N 

We note here that as N — > 00, t]n — > 1 and rjv — > 00. 
Both results are very unusual, since typically 77-exponents 
are small. For reference, in the O(N) model 77 — > in 
the N -> 00 limit and r\ w 0.037 for the 0(N = 3) 
model |l62l |. 

We now turn back to the study of the lattice designer 
J1-J2-Q Hamiltonian. By analyzing the size dependence 
of the correlation functions in Eq. (|23p at the location of 
the critical points shown in Fig. [51 it is possible to esti- 
mate values for 77 at and r\y as a function of N [13, Il37l |. 
Results for 2 < N < 12 are summarized in Fig. [7J The 
lattice Hamiltonians reproduce the "smoking gun" DQC 
feature that 77 at — > 1 and rjv — > 00 in the N — > 00 limit. 
Furthermore, extrapolations of the data, with the con- 
stant terms subtracted off, show quantitative agreement, 
within a few percent, with the corrections in the large- TV 
expansions, Eq. (|2~4"]) . Note that the correct description 
for large N is far from trivial, since the continuum the- 
ory was not derived from the microscopic Hamiltonians 
considered here but was posited as a generic theory of 
the Neel VBS transition. 



The fact that the exponents agree so well lends strong 
positive support in favor of the DQC scenario for large 
N . Clearly, it would be very interesting to compute fur- 
ther 1/N corrections analytically, to test the theory in 
this quantitative way down to the smallest JV. We note 
that the exponent 77^ for the SU(2) case has been esti- 
mated based on classical Monte Carlo simulations of the 
non-compact CP 1 model, and also of the 3D Hciscnberg 
model with "hedgehog" topological defects supp ressed 
(which should lead to the same universality) |l0l| . The 
exponents are in rough agreement with those of the SU(2) 
J-Q model, but it would be very interesting to push these 
classical simulations to higher precision as well. Another 
study of a classical model argued to realize the DQ C ac- 
tion was found to have a first-order transition 11631. 



E. SU(A0 Bilayer 

The study of the destruction of SU(JV) magnetic order 
in the bilayer geometry, illustrated in Fig.|SJa), provides 
an interesting example to test our understanding of the 
role of Berry phases at DQC points, since these should 
cancel between the two square lattice layers. Consider a 
bilayer in which each layer is described by some combi- 
nation of the J1-J2 interactions discussed in the previous 
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section and with the following coupling between the lay- 
ers, 



SPIN LIQUIDS AND DECONFINEMENT IN 
U(l) MODELS 



Hj 



(25) 



where the sum over [ij] is taken between spins vertically 
next to each other in the bilayer (i.e., [ij] form dimers, 
as in Sec. IIII A[) . For any N, the J^-only model has a 
featureless non-degenerate ground state consisting of a 
product of vertical valence bonds (singlets) — the SU(A) 
generalization of the VBL state discussed in Sec. IIII Al 
The introduction of small intra-layer couplings like the 
J\-J<i interactions will make the ground state deviate 
from a simple product state by increasing the amount 
of entanglement, but these couplings cannot destabilize 
the VBL paramagnet, as long as they are small compared 
to J x ; see Fig.^b). 



1. Neel-VBL transition 

When the value of J2/J1 places the system within the 
Neel phase when J± = 0, there is a transition from the 
Neel state into the VBL for increasing J±. This transi- 
tion in the SU(2) case is continuous and belongs to the 
3D 0(3) universality class, as d iscussed in Sec. IIII Al For 
N > 3, QMC calculations jl64l | show that the transition 
is first-order, however, as also predicted by a mean-field 
theory. For N = 3 the transition appears continuous 
based on the largest system sizes studied. If this is in- 
deed the case, it would be in the universality class of the 
compact CP 2 model j4f|. Further studies of the N = 3 
model should be carried out to resolve this issue. 



2. Neel-VBS transition 

A very interesting aspect of the SU(iV) bilayer mod- 
els is the limit where J± is weak and there is a Neel- 
VBS transition as a function of J2I J\- It is well known 
that the Berry phases cancel in the long- wavelength limit 
between the layers. Nevertheless, both Neel and VBS 
phases must clearly be stable to a small finite interlayer 
coupling. What is the fate of the J± = deconfined 
Neel-VBS transition in the presence of the bilayer cou- 
pling? Remarkably, the QMC simulations show that the 
transition between the very same Neel and VBS phases 
becomes first-order in the bilayer geometry. This find- 
ing can be understood as a restoration of the Landau 
paradigm, due to the cancellation of Berry phases which 
is a r elevant perturbation at the deconfined critical point 
[l64| . This is another piece of strong evidence for the 
correctness of the DQC theory, as well as the more gen- 
eral field-theoretical understanding of the role of Berry 
phases. 



One of the main foci of research into designer Hamil- 
tonians is the search, detection and characterization of 
quantum spin liquid (QSL) phases. That is, although 
there has been much progress in understanding prop- 
erties of QSL states from effective field theories, there 
are very few microscopic models amenable to unbiased 
calculations that can be argued to exhibit a QSL with- 
out controversy. The recent spate of high-profile numer- 
ical stu dies i dent ifying possible QSL candidate Hamil- 
tonians |l65l - ll67| has sparked intense activities on the 
theory side, highlighting the synergy between analytical 
theory and numerics. Thus, we are interested in finding 
the simplest models harboring QSL states, that are also 
amenable to large-scale simulation by QMC. 

Much of the reason for the scarcity of microscopic mod- 
els is that fundamental ingredients required to give rise 
to QSL phases, namely geometrical frustration, are also 
typical harbingers of the sign problem. Despite this, the 
absence of the sign problem docs not fundamentally pre- 
clude the existence of spin liquid physics. This is most ev- 
ident in Hamiltonians with U(l) symmetry — anisotropic 
spin (or Bose-Hubbard) models specifically designed to 
capture the essential ingredients needed to promote a 
QSL phase without causing the sign-problem. A cru- 
cial fact in this regard is that there are no limitations on 
the diagonal couplings, because the sign problem is only 
caused by the off-diagonal matrix elements in the chosen 
QMC basis (see Sec. [TT|) . Thus, diagonal frustration can 
be simulated in any form. 

Examples of both 2D gapped [16^171 1 and 3D and 
gapless [l72| QSLs have recently been studied with large- 
scale QMC simulations of such sign-problem free U(l) 
models. In this section, we will explore one class of 
Hamiltonians that contains the simplest type of QSL 
phase, which can be related to a Z2 gauge theory We will 
demonstrate that universal topological properties calcu- 
lated in QMC can be compared to this effective field- 
theoretic description. In addition, an associated critical 
point separating an XY-ferromagnet and a QSL can be 
demonstrated to be of an exotic "fractionalized" type, 
based on comparisons of universal quantities calculated 
from QMC and a related quantum field theory 



A. BFG Hamiltonians with Z2 spin liquid phases 

One successful recipe for designing Hamiltonians with 
QSL phases amenable to QMC study was pioneered by 
Isakov and co-worker s for certain B ose-Hubbard models 
on the kagome lattice jl68l . ll69lll7ll |. Models in this class, 
first proposed analytically by Balents, Fisher and Girvin 
[l73j (BFG), are constructed partly based on their rela- 
tion to quantum dimer models (QDMs), where a strong 
local constraint (the number of dimers emanating from 
a site being fixed) is coupled with a quantum tunneling 
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FIG. 9. (a) Spins (black dots are S z = 1/2, empty sites 
S~ = —1/2) in a groundstate configuration satisfying Hq, 
Eq. ()26[) . on the kagome lattice. The degenerate manifold of 
such configurations maps onto a triangular-lattice (dashed) 
classical three-dimer configuration, wh ere e ach spin-up corre- 
sponds to a dimer on the dual lattice |l73l ]. (b) A plaquette- 
flip -ff r ing operation preserves the cluster-charging constraint, 
and can be seen to correspond to a rhombus-flip of dimers 
on the dual lattice. A defect is created by flipping one spin, 
marked "+" in (a). This defect can be thought of as associ- 
ating with two hexagonal plaquettes, indicated in grey in (b), 
each with two spin-up per hexagon, corresponding to decon- 
fined fractionalized spinons. 

that allows the state to fluctuate, while keeping the con- 
straint satisfied. QDM Hamiltonians can be constructed 
where the groundstate superposition of dimer configu- 
rations breaks no symmetry — being a type of QSL (or 
resonating valence bond) phase, that has an eme rgent 
gauge symmetry and related topological order |l74l |. 

The recipe to construct QSL phases in XXZ spin (or 
boson) models is to map the dimer constraint to a plaque- 
tte spin configuration on the dual lattice; the dimer may 
represent an S z spin configuration on a bond which is 
frustrated (or unsatisfied) . Then, the role of the quantum 
dimer tunneling moves is done by a spin-exchange pro- 
cesses (which may not be a simple two-body term). Such 
models have been called cluster- charging models |169| . 
because spin configurations summed over a local clus- 
ter (e.g. a lattice plaquette) are penalized when differing 
from some chosen value. For example, on the kagome 
lattice (Fig.EJ), 

H = VY / (S& 2 , (26) 

o 

where 

So = E s t ( 2? ) 

is a cluster-charging potential term which favors three 
spin up and three spins down per hexagonal plaquette. 
This Hamiltonian alone will promote a highly degener- 
ate manifold of groundstate configurations, where each 
kagome-lattice hexagon satisfies the plaquette constraint 
but is otherwise disordered. Like all operators which are 
entirely diagonal in the chosen QMC basis, H does not 
affect the simulation with a sign problem, no matter how 



the interaction is frustrated, making it a suitable ingre- 
dient for a designer Hamiltonian. 

In order to promote a QSL phase from this classical 
degenerate manifold, one desires to add quantum fluctu- 
ations to the Hamiltonian such that the cluster-charging 
constraint is not violated. In this case, a ring-exchange 
term, 

ffri„ g - -K ^(S+SrS+Sf + h.c), (28) 

which operates around a "bow-tie" plaquette on the 
kagome lattice (Fig. promotes quantum fluctua- 
tions while preserving the cluster-charging constraint 
Hq. Thus, one expects the kagome lattice model with 
H = Hq + H ring to retain a disordered groundstate with 
quantum fluctuations — a good candidate for a quantum 
spin liquid phase. In addition, with K > , such terms 
do not have the sign problem for QMC [l75j . 

BFG [173| | were the first to show that variations of this 
kagome-lattice model support the simplest type of Z 2 
spin liquid — in part icular, at an exactly soluble Roksar- 
Kivelson point jl3lj . Subsequently, using QMC simula- 
tions, Isakov and co-workers showed that a variant of the 
cluster-charging Hamiltonian with different short-range 
spin exchanges, 

Hl, 3 = -il, 3 Yl + ( 29 ) 

(y)i,3 

support a Z2 spin liquid phase. This exchange term can 
either connect the first, second an d third neighbors (ij) 3 
on the kagome lattic e |l68l . Il69j |; or it can be nearest- 
neighbor only, (ij)i [l7l| . Again, for ti t s > no sign 
problem exists. In addition to H = Hq + H3 and H = 
Hq + H\, the Z2 spin liquid has been demonstrated on 
H = Hi + -f/ring, where the cluster-charging term is not 
explicitly required and the Ham ilton ian consists of two 
competing kine tic-e nergy terms |l70| . sometimes called 
the J-K model p75| . 

In the next section, we explore one procedure by which 
QMC simulations can positively identify such spin liquid 
phases. In section IV C[ we examine the XY-ferromagnet 
(or "supcrfluid" in the boson language) to QSL phase 
transition, which is common to these models, as a dif- 
ferent example of a deconfined or fractionalized quantum 
critical point. 



B. Identifying the Z2 spin liquid with entanglement 

The smoking gun signature for a gapped spin liquid 
phase could be argued to be the identification of an emer- 
gent gauge theory, and the associated topological order 
[l76j |. One manifestation is the existence of a ground 
state degeneracy, which is topological in nature, in the 
simplest (Z2) case being four-fold on a torus. Previous 
QMC studies of QSL states have demonstrated the dif- 
ficulty in clearly identifying this degeneracy, due to the 
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tendency for tunneli ng b etween the equal-energy states 
on finite-size lattices [l68l | . Similarly, recent DMRG work 
studying candidate gapped spin liquid states have been 
unable to identify the expected topological degeneracy 

[III Gil. 



1. The boundary law 

Luckily, there is another tool suited for the positive 
identification of this topological order, through measure- 
ments of the Renyi entanglement entropy, introduced in 
section Hi Gl The quantities S a are well studied in quan- 
tum information science; they quantify correlations for 
a system in a basis-independent way, which makes them 
suitable in particular for characterizing phases that do 
not have an explicit broken symmetry (i.e., QSL phases). 
In fact, some authors have championed the entanglement 
entropy as a paradigmatic analog to symmetry breaking 
for non-Landau phases; here, the concept of Long Range 
Entanglement (LR E) replaces the concept of Long Range 
Order (LRO) [1761 ]. For example, in a gapped phase with 
LRO (such as the VBS phases mentioned above), the en- 
tanglement entropy is expected to obey the "boundary 
law" (also called the "area law" ) 



S a — 



(30) 



where a is a non- universal constant, t is the length of 
the boundary between regions A and B, and terms de- 
noted by • • • scale away at least as fast as 0(l/£). This 
can heuristically be considered as "short-range" entan- 
glement. In contrast, in a gapped spin liquid phase (with 
no broken symmetries) scaling of entanglement entropy 
is predicted to be of the form 



S a = at - 7 + 



(31) 



Here, the subleading term to the boundary law is a uni- 
versal conection called the topological entanglement en- 
tropy [l77l Tl80j. In a gapped quantum spin liquid phase, 
it is related to the emerge nt g auge symmetry, indepen- 
dent of the Renyi index a [18l[ . In p artic ular, in the Z2 
spin liquid state, 7 = ln(2) at T = 179]. The addition 
of this universal subleading constant to the boundary law 
motivates the notion of LRE for a gapped quantum spin 
liquid phase. 



2. Calculating the topological entanglement entropy 



Using the QMC estimators described in Section III Gl 
it is possible to isolate the topological entanglement en- 
tropy from the leading-order boundary law, as well as 
subleading corrections due to effects such as corners, by 
taking measurements using several different geometries 
A and performing additions or subtracti ons which iso- 
late only 7. For example, the Levin-Wen [179| construc- 
tion isolates 27 from Eq. (|31[) . requiring measurement of 



four unique region A geometries. For a BFG Hamilto- 
nian discussed in the previous section, it is possible to 
extract the topological entanglement entropy in a QMC 
simulation using t his c onstruction. As noted by Castel- 
novo and Chamon [l82| through direct calculation on the 
toric code (which realizes a Z2 fractionalized spin liq- 
uid groundstate) , the value of 27 is approached through 
two temperature crossovers on a finite-size system, re- 
lated to the energy of the quasiparticle excitations in the 
groundstate — a Z2 charge corresponding to the spinon, 
and a Z2 flux corresponding to the vison. Each excitation 
contributes a value of ln(2) to the topological entangle- 
ment entropy. As illustrated in Fig.[T0l the realization of 
one of these plateaus at finite-temperature gives positive 
identification that the groundstate of the model is indeed 
in a topological Z2 QSL phase. 



C. Deconfinement at the XY* transition 

The fractional spinons of the Z2 spin liquid phase (il- 
lustrated in Fig. O can be thought of as a "square-root" 
operator of the physical spins. Employing the mapping 
from S = 1/2 spins to hard-core bosons, we can regard 
the physical spin operator S + or (equivalently) boson 
operator = cj^cj^ as being composed of two creation 
operators 4> for the fractionalized spinon particles. 



1. The anomalous dimension 

Surprisingly, it is predicted that the quantum phase 
transition between the XY-ferromagnet and the Z2 spin 
liquid discussed in the last section can actually be me- 
diated not by the physical bosons (which would undergo 
a quantum phase transition in the 3D XY unive rsality 
class) , but by the fractional spinon fields dH-EHl. This 
leads to the same exponents z and v as in the ordinary 
XY critical point. However, the exponent 77 controlling 
the equal time correlation function (6t(0)&(x)) is mod- 
ified. In a way similar to the DQC point discussed in 
Sec. IIV C[ 77 = 1 in the TV — > 1 limit and is expected to 
be anomalously large also for TV = 1. 

Previous estimates for 77 were obtained from field- 
theoretic treatments. From an 1/TV expansion, Ref. jl84l ] 
obtains 



77=1 



32 
3^TV' 



(32) 



which gives 77 = 2.08 for TV = 1. A more accurate 
value was obtained from a combination of field theory 
and Monte Carlo simulations of the co rrelation of a com- 
posite operator in the 3D XY model [l86l Il87| , leading 
to 7/ w 1.47(3). 

The measurement of 7; in the kagomc BFG model 
is challenging, requiring measurement the equal time 
Green's function in real space: G(x) = G(r ~ 0, r) = 
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FIG. 10. Main graph: Temperature dependence of the topo- 
logical entanglement entropy in the Z2 QSL phase (V/ti = 8) 
of the Hamiltoni an H p + Hi for three different L x L x 3 
kagome lattices |l7l| . The quantization of the plateau at 
ln(2) (dashed line) identifies the emergent gauge symmetry; 
one expect another temperature crossover to 27 = 21n(2) at 
some lower temperature. Inset: finite-size scaling of G(L/6) 
at the XY* quantum critical point (V/ti = 7.0665), which 
scales as 1/L 2 ' 9 3 (st raight-line fit) giving 77 = 1.493. Repro- 
duced from Ref. [I89f ] . 

(&t(0)6(r)) [H GUI, done by keeping track of the de- 
fects created in the SSE directed-loop algorithm [67} as 
it traverses the (d+l)-dimensional space-time QMC sim- 
ulation cell. At the critical point, G(r) should decay as 
l/r 1+n , and finite size effects can be minimized, e.g., by 
measuring G(L/6), which decays as l/L 1+T >, as a function 
of L. Fig.fTUlshow an algebraic decay with rj = 1.493(10); 
a value consistent with 77 for the composite operator in 
the 3D XY field theory mentioned above. The excellent 
agreement between the QMC simulation results and the 
composite-boson field theory is a remarkable confirma- 
tion of fractionalized universality in this model; syner- 
getic work between QMC and quantum field theory is 
uniquely capable of positive identification of the exotic 
quantum phase transitions in this model. 



2. Detecting fractionalization 

As discussed in Ref. |l89j |. there are consequences of 
this novel fractionalization in other universal quantities. 
Remarkably, even though the fractional spinons become 
gapless at the quantum critical point, there are still ob- 
servable remnants of the topological nature of the Zi 
spin liquid. This can be seen by looking at the winding 
number distribution at the XY* point, and comparing 
it to the winding number distribution at a regular 3D 
XY transition. Namely, since the fractional spinons are 
bosonic, they can be modeled by a "simulation within a 
simulation" , achieved by projecting physical bosons at an 
XY transition to wind strictly in pairs around the space- 
and time-periodic simulation cell. Since this distribution 



is universal, it is expected to be the same (up to a non- 
universal velocity) as the winding number distribution 
of the physical (composite) boso ns in the XY* transi- 
tion. Observation of this fact [l89j | is a striking confirma- 
tion that topological properties can still be present even 
when fractional excitations are gapless. These universal 
winding number distributions are possible to compute in 
principle in field theory; hence they are another example 
of the synergetic connection between and quantum field 
theory and QMC simulations on designer Hamiltonians. 

VI. DISCUSSION 

We have presented several examples of the use of 
judiciously constructed designer Hamiltonians for test- 
ing field theories and exploring various quantum many- 
body phenomena related to quantum-criticality and ex- 
otic ground states. Here we discuss some further generic 
aspects of this approach, followed by a brief discussion of 
other interesting related studies of systems not covered 
in this Review, and an outlook on future directions. 

A. Designer Hamiltonian, stoquasticity, and the 
sign problem 

The concept of designer Hamiltonians that we have 
introduced he re is different from that of "stoquastic" 
Hamiltonians jl9C| . which are those for which all off- 
diagonal matrix elements are non-positive (i.e., sign- 
problem free) in the standard basis. This is normally the 
^-component basis for spin systems, although a broader 
definition with respect to an arbitrary basis has also been 
presented [191], which would seem to indicate that any 
Hamiltonian could in principle be cast into the stoquastic 
class. Our criterion is more one of scientific utility than 
a strict mathematical definition — models that are free 
from sign problems cither automatically (in a practically 
useful basis), or through some other way of circumvent- 
ing it (e.g., the Meron algorithm [5Fj| or determinant- 
bas ed QMC of particle-hole symmetric Hubbard models 
[49L I192L Il93| ) . To be considered a designer Hamiltonian, 
such a model should also represent a prototypical case of 
some interesting quantum-man body phenomenon which 
is difficult to access in an unbiased manner in other ways, 
i.e., the model is designed to shed light on this prob- 
lem (and must be "de-signed" to be practically useful for 
studies on large scale). 

1. Limitations of designer Hamiltonians 

An important fundamental problem is whether it is al- 
ways possible to construct such designer Hamiltonians, or 
whether certain types of states are beyond reach in prac- 
tice because their physical properties arc fundamentally 
tied to the difficulties in circumventing the sign problem. 
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For instance, thc"Bose metal" spin- liquid phase with a 
Fermi sea of emergent fermions appears to have a com- 
plex sign structure that would be difficult to repr oduce 
in a sign- free Hamiltonian in a simple basis [194j . On 
the other hand, as we have reviewed in Sec. IV Al Z2 spin- 
liquid ground states have been found in models with U(l) 
symmetry. An interesting open challenge is to fin d the Z2 
spin liquid (which itself is a broad class if states |l95| ) in 
sign-free SU(2) invariant models. It is certainly possible 
that these states cannot be realized with sign-problem 
free Hamiltonians in a simple basis, but one can still not 
exclude that designer Hamiltonians can be constructed 
for which some useful basis can be found, or for which 
the sign problem can be solved in some other way. In- 
deed, the recent discovery of a spin liqu id in the Hubbard 
model on the honeycomb lattice 1661 ] seems to provide 
an example, although it remains to be seen whether this 
state really is a Z2 topological state. 



2. Computational complexity 



Re gard ing the sign problem itself more broadly, 
Ref. |l96| presented a proof that the sign problem can- 
not be solved in general. The problem here was defined 
in terms of constructing a QMC algorithm for which the 
computational effort scales as a power-law in the system 
size. A particular problem is considered for which the 
removal of the negative signs still leaves a problem for 
which the Monte Carlo sampling is believed to scale ex- 
ponentially in the system size and, thus, the sign prob- 
lem was not solved, only recast into a different form. 
However, the model system chosen in this demonstra- 
tion was that of a spin glass, and the conclusion reached 
is equivalent to the statement that the quantum spin 
glass is no easier to solve than the corresponding clas- 
sical spin glass. This is hardly a surprising result. The 
proof does not address QMC solutions of systems that 
are not directly associated with the well known difficul- 
ties of simulating classical frustrated "glassy" systems 
with complicated energy landscapes. Whether or not the 
sign problem can be circumvented (at least for represen- 
tative designer Hamiltonians) for translationally invari- 
ant systems such as fermionic Hubbard models without 
particle-hole symmetry and SU(2) frustrated quantum 
spins models remains an important open question. 

Some non-trivial solutions of sign problems are worth 
mentioning. In some cases, local basis transformations 
can be used to r ender all off diagonal matrix elements 
non-positive |197( . One can some times apply more so- 
phisticated ways to circumvent the problem, e.g., with 
the Meron alg orith m [55j j or the newly discovered fermion 
bag approach jl98l |. It remains to be seen whether these 
exiting developments can be applied to a wider range of 
fermionic and frustrated designer Hamiltonians. 



B. Other systems and future prospects 

1. Spin models 

The models and quantum states we have discussed here 
only represent a small sample of recent research into fas- 
cinating topics in quantum magnetism. Other notewor- 
thy examples include syst ems exhibiting Bose-Einstein 
condensation of magnons jl99j |. spin supersolids j20C| . 
qud rupo lar order [20 1[ |202| , and quantum compass mod- 
els |203| |. Here as well, QMC studies of designer Hamil- 
tonians played a crucial role in uncovering novel physics, 
and unresolved issues will require the kind of synergistic 
interplay between low-energy theories and QMC simula- 
tions that we have emphasized in this Review. 

Within the classes of spin models we have discussed 
here, although the essential physics now has been ascer- 
tained, there still remains important work to be done to 
improve on the quantitative aspects, e.g., the critical ex- 
ponents. In the case of the studies aiming to connect 
SU(AT) QMC simulations with large- N expansions, sum- 
marized in Sec. lIVDl and Fig. [71 it would be very useful to 
go to still larger lattices, to check the convergence of the 
computed exponents more precisely. This is important, 
in particular, in light of the fact there there are significant 
scaling corrections in quantities such as the spin stiffness 
for N — 2. Such corrections could potentially also affect 
the fitting procedures underlying the exponents graphed 
in Fig.[JJfor small N. When higher-order 1/N analytical 
results eventually become available, it will be important 
to know the values of the exponents without possible re- 
maining effects of scaling corrections. What is needed 
here is high-accuracy (small error bars) results for the 
correlation functions on large lattices. It will also be 
very interesting to use the methods discussed in Sec. Ill Gl 
to quantify the degree of entanglement in the SU(2) and 
SU(A0 models. 



2. Fermions 

Another very interesting aspect of the designer Hamil- 
tonians we have discussed here is that they can also be 
doped with fermions. While this introduces a sign prob- 
lem, it should be possible to study large lattices doped 
with a small number of holes, since the magnitude of the 
sign problem only depends on the number of fermions. 
This low-density limit of the SU(iV) models should al- 
ready be sufficient for making connections with large- 
theories of "holonic" or "strange" metallic states, that 
have be en pr edicted to arise out of the DQC system upon 
doping [l60j |. 

An interesting conceptual question is whether one can 
study bosonic systems in which fermions are emergent 
particles. A positive answer may be suggested by the 
demonstrated existence of the Z2 fractional phase in a 
purely bosonic model, as we discussed here in Sec. IV Bl 
The finite-temperature structure of the topological en- 
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tanglcment entropy identifies that both a charge (a 
spinon) and a Z-i fl ux (a vison) are realized as excitations 
in this model |l7l| . This raises the possibility that, with 
a suitable "binding" potential {204} coupling these two, 
an extension of the model with an emergent fermionic 
excitation might be created fl76j . The ability for QMC 
to measure the entanglement entropy in a state with an 
emergent fcrmion could also have wide implications for 
the purported application of tensor network algorithms 
as a "solution" to the sign problem, if such a state were 
found to violate the boundary law, Eq. (j30l) . 



3. Disordered systems and impurities 



can observe a vortex in the ord er param eter, and at the 
SU(2) and SU(3) critical points [l3l[l53, where impuri- 
ties provide an interesting way to probe the anomalously 
large corrections to scaling in these systems. Impurities 
in ID chains (with open edges) have also been studied 

In general, impurities and disorder in interacting quan- 
tum systems present rich opportunities for further ex- 
ploring exotic states and phase transitions. Synergistic 
combinatio ns of low-energy fiel d theories. RG schemes for 
impurities [223| and disorder (2l l [223 - 1226} . and QMC 
studies of designer Hamiltonians will be needed to make 
progress here. 



In this Review we have focused on translationally in- 
variant systems. Systems with disorder, i.e., random 
couplings, dilution, etc., is another area where designer 
Hamiltonians can provide unique insights. To men- 
tion one example of a spin system for which the stan- 
dard classical-quantum correspondence fails, in the 2D 
5 = 1/2 Heisenberg antiferromagnet diluted with non- 
magnetic impurities, the interplay of classical percola- 
tion an d qu antum fluctuat ions lead to low-energy exci- 
tations [l47{ not captured 205 within the conventional 
d — > (d + l)-dimcnsional mapping. Other unanticipated 
effects of dilution on quantum-critical scaling have also 
been obser ved that s till lack explanation and deserve fur- 
ther study 206l4210j [. For the designer Hamiltonians dis- 
cussed in this Review, disorder effects have not yet been 
studied. One can anticipate interesting behaviors in VBS 
states, at the Neel-VBS transition, and in the Z2 spin liq- 
uids. 

Beyond disordered spin models, various manifestations 
of the " dirty boso n" problem remain intriguing and chal- 
lenging [21 lM213 1 . and provide excellent examples of the 
indisp ensa ble role of QMC studies of designer Hamilto- 
nians [214j . An example is the extended hardcore Bose- 
Hubbard model with random frustrated interactions on 
a 3D cubic lattice. QMC simulations have been instru- 
mental in discovering that glassiness can coexist with su- 
perfluidity, and that quantum fluctuations and random 
frustration are both cruc ial ingred ients needed to stabi- 
lize a "superglass" state [2151 . 1216} . 

In addition to extensive disorder, single-impurity prob- 
lems are also interesting and important. Effects of a sin- 
gle missing or added spin have been investigat ed with 
QMC in the 2D and 3D Neel states [217^2191 and at 
con vent ional quantum-critical points in dimerized mod- 
els j220} . QMC studies have also been d one o n vacancies 
in the VBS state of the 2D J-Q model [22l} . where one 



4- N on- equilibrium dynamic scaling 

A recent development on the methods front is the 
adaptation of QMC algorithms such as those discussed in 
Sec. [TT] to non-equilibrium setups [Ti} . By studying time 
evolution out of a non-equilibrium state in imaginary 
time, standard sign-problem free systems can be analyzed 
in novel ways. For instance, it was found in Ref. fl4} 
that the dynamic scaling properties when quenching to a 
quantum-critical point is the same in real and imaginary 
time, and that such protocols are useful for extracting 
the dynamic critical exponent. Moreover, the fidelity and 
full geometric tensor, which characterize the geometrical 
properties of the state space, can be accessed this way. 
One can anticipate many fruitful applications of this un- 
conventional QMC approach, which gives access to dy- 
namics beyond the commonly used numerical analytic 
continuation of QMC-co mput ed equilibrium imaginary- 
timecorrelation functions |227} . 
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